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Abstract 

The discovery of neutrino masses has provided strong hints in favor of the possibihty 
that B-L symmetry is an intimate feature of physics beyond the standard model. I 
discuss how important information about this symmetry as weU as other scenarios for 
TeV scale new physics can be obtained from the baryon number violating process, 
n — n oscillation. This article presents an overview of different aspects of neutron-anti- 
neutron oscillation and is divided into the following parts : (i) the phenomenon; (ii) 
the physics, (iii) plausible models and (iv) applications to cosmology. In particular, it 
is argued how the discovery of n — n oscillation can significantly affect our thinking 
about simple grand unified theory paradigms for physics beyond the standard model, 
elucidate the nature of forces behind neutrino mass and provide a new microphysical 
view of the origin of matter in the universe. 



1 Introduction 



Particle oscillations are familiar phenomena in both classical and quantum mechanics and 
have provided a wealth of information about the nature of matter and forces acting on them. 
The most elementary example of a pendulum oscillation is caused by the force of gravity 
and therefore provides information about the strength of gravity. In the domain of quantum 
mechanics, two states close by in energy can transmute into one another if the Hamiltonian 
for the system includes a force that connects them and if certain coherence conditions are 
satisfied. Turning this question around, one can conclude that if experimentally an oscillation 
between states, forbidden in the standard model is observed, it will then imply the existence 
of hitherto unknown classes of interactions. Many such instances in the domain of particle 
physics are known where study of particle oscillations have provided important clarification 
about the fundamental nature of underlying physics. Kaon and neutrino oscillations are but 
two glowing examples of this. The first in combination with the recently observed B — B 
oscillations confirmed the nature of CP violation in standard model. The latter has for 
the first time revealed the existence of new physics as well as new symmetries beyond the 
standard model. 

The question then arises as to whether there are other such possibilities that need to be 
experimentally explored to probe physics beyond the standard model (SM). In this review, 
I argue that neutron-anti-neutron (n — n) oscillation provides one such unique example and 
is extremely timely in the aftermath of one of the fundamental discoveries in our field i.e. 
the neutrinos are massive and can provide important clues to our understanding of new 
symmetries and forces behind the neutrino mass generation. 

The n — n oscillation[l] is a unique kind of oscillation phenomenon compared to kaon 
oscillations in that it breaks one of the once sacred conservations laws of physics i.e. con- 
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servation of baryon number which keeps all matter stable (just like kaon oscillation breaks 
strangeness quantum number). It also effectively makes neutron a Majorana fermion (albeit 
with a very tiny Majorana component), a point which is of some historical importance since 
in the original paper of Majorana, he contemplated neutron being the particle which is its 
own anti-particle. 

The fact that n — n oscillation breaks baryon number requires us to first ascertain that 
the vast wealth of information available about the stability of matter from experiments and 
general cosmological observations allow for this process to have a strength that will make it 
observable with current technology. We show below that this is indeed the case. 

Secondly as far as the physics implications go, another class of processes that probe 
baryon number nonconservation is proton decay [2, 3]. They have been the focus of exper- 
iment as well as theory for two reasons: first of all, simple grand unified theories [4] based 
on SU(5) and SO(IO) predict its existence at a level not far from experimental capabilities; 
secondly, it used to be thought in the early 80's that these theories may also explain the 
origin of matter using proton decay as one of its key ingredients. 

This situation however changed drastically with time due to two developements. 

(a) The discovery of neutrino masses was confirmed in 1998 and a popular paradigm for 
the origin of their smallness became the seesaw mechanism[5] , where new symmetry, B-L 
respected by then standard model forces had to be broken^ . Note that the dominant proton 
decay modes of interest p e^Ti° respect the B — L symmetry and is therefore nor suited 
for the study of forces responsible for neutrino mases. On the other hand, n — n oscillations 
break B-L symmetry by two unjits exactly like the seesaw mechanism for neutrino masses 
does and therefore is uniquely suited for the study of the detailed nature as well as the 

^Without any new symmetry, the seesaw scale would naturally be expected to be the Planck scale, which 
yields too small neutrino masses. A simple way to understand why the seesaw scale is lower than the Planck 
scale is to assume that it is associated with a symmetry e.g. B-L . 
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scale of seesaw mechanism. There is a large class of interesting gauge models where indeed 
neutrino Majorana masses via the seesaw mechanism leads directly to n — n oscillations. 

(b) The second development is that a new mechanism for understanding the origin of 
matter was proposed in mid 80's that did not use proton decay but rather used the seesaw 
mechanism for neutrino masses[6] to generate matter-anti-matter asymmetry where again 
B-L symmetry played a key role. Although grand unified theories such as SO(IO) could 
incorporate the B-L symmetry, these two developments took a lot of the wind from the 
argument in favor of proton decay being necessarily the most theoretically relevant mode for 
baryon nonconservation. 

Furthermore, when the seesaw mechanism is embedded into GUT theories that include 
B-L symmetry, they generally tend to fix the scale of B-L symmetry breaking to be at the 
GUT scale, while there is no compeUing reason for the seesaw scale to be that high. It could 
even be at the TeV scale without making the theory too unnatural. The seesaw scale could 
not therefore be probed by searching for proton decay whereas as already noted, in a large 
class of interesting theories, n — n oscillation does. 

The question then arises as to whether all physics beyond the standard model is encap- 
sulated in grand unified theories with neutrino masses being also a signal of the same GUT 
scale or whether the neutrino masses indicate an alternate path that contains rich physics 
far below the GUT scale. As we show below, n — n oscillation provides an effective probe 
of the second possibility which envisions a rich new layer of physics at intermediate scales 
that involves low B-L breaking. Theoretical models indicate that a wide range of scales, 
anywhere from a few TeV to 10^^ GeV as well as the associated new physics could be probed 
by pushing the n — n oscillation oscillation time by three orders of magnitude over the cur- 
rent experimental lower limit. In particular, its discovery will rule out simple grand unified 
theories or in the very least have significant affect on their detailed nature. Furthermore 
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it will have profound effect on the nature of physics beyond the standard model as well as 
provide a new microphysics view for the origin of matter, different from leptogenesis. It also 
provides rich new physics that can be probed at the Large Hadron Collider. 

The bottom line is therefore that it is urgent to pursue a dedicated search for for n — n 
oscillation that can push the sub-GUT physics frontier to a new level. 

This review is organized into the following parts: in the first part, I discuss the phe- 
nomenon, its connection to nuclear instability ; in the next section, we discuss the broad 
character of microphysics that would lead to n — n oscillation and the kind of scales that can 
be probed by this. In the next section, the nature of physics beyond standard model probed 
hj n — n oscillation will be discussed including a discussion as to why B-L is an important 
symmetry to be probed. In this section, we discuss specific and plausible gauge models for 
n — n oscillation both with and without supersymmetry. In the final section, we discuss the 
implication of observable n — n oscillation for baryogenesis and show how one can have a 
new class of models for late baryogenesis after electroweak phase transition and sphaleron 
freeze-out. 



2 Phenomenology of n — n oscillation 

In order to discuss the main features of n — n oscillation in various environments, it is 
sufficient to consider a simple 2x2 Hamiltonian for the evolution of an initial beam of slow 
moving neutrons and evaluate the probability of finding an anti-neutron in this beam after 
a time t. 

= ( En Sm\ fn\ 

where 5m denotes the n — n mixing which parameterizes at the nucleon level the underlying 
physics that breaks baryon number by two units. We will give examples of models where 
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this can happen at the appropriate level in a subsequent section. Note that we have left the 
energies of the n and n arbitray since they could be affected in different ways in the presence 
of an external environment such as the nuclear field or a magnetic field etc. 

The exact expression for the probability to find an antineutron at time given that at 
t = 0, P„(0) = 1 and P^(0) = is: 



45m 
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where = — En. The fact that neutron decays does not have any effect on this 

formula as long as the neutron fiight time is much smaller than the life time of neutron. 

Let us consider two extreme cases of this formula which also happen to be interesting for 
experimental purposes. For Sm <^ AE, this expression reduces to: 

Pn ~ sin'AEj (3) 

There are two special cases of this realizable in nature: 

• Case (i): AEt <^ 1: In this case we have, Pn-*n ~ ((^"7. • = This case 
corresponds to free neutron oscillation in vacuum. 

• Case (ii): AE • i ^ 1: Pn^n ~ h ( oaf ) T^^is for instance corresponds to bound 



2 \2AE„ 

neutrons inside nucleus "oscillating" to anti-neutrons . 



2.1 Stability of Nuclei and limit on r, 



n—w 



As with all baryon number violating interactions, the existence oin — n oscillation will lead 
to nuclear instability. This will happen via a conversion of neutron to an anti-neutron inside 
the nucleus which will then annihilate with the surrounding nucleons and lead to instability. 
Clearly, if the nuclear transition probability inside the nucleus and in vacuum were equal. 
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this would cause all matter to disappear in a much shorter time than the age of the Universe 
and would be unacceptable. However, it is an experimental fact that inside the nucleus the 
neutron and the anti-neutron experience different nuclear potentials. In fact the difference 
is so large {E^ — Ef^ ~ 100 MeV or so) that the transition probabihty formula relevant for 
this case is that in case (ii) above and leads indeed to nuclear instability lifetimes of order 
10^^ years or so. Careful and detailed analysis of this question has been done by Dover, Gal 
and Richards[ll], Alberico et al[12] and more recently by Gal and collaborators [13]. Below 
we present a crude representation of the above works to get the basic physics across: 



Pinuc. o 



/ 6m 
[2AK 



(4) 



where Pn^n represents the free neutron oscillation probability in vacuum and 5m is the 
AB = 2 off diagonal matrix element in the n — n mass matrix in Eq. (1) and is related 
to the oscillation time Tn-n — 2iTkn- P^^sent ILL limit on T„_fj > 10^ sec. translates to 
Sm < 6 X lO'^^ GeV. A crude estimate of the nuclear instability hfe time in the presence 
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n — n oscillation from these intuitive considerations is TnucI. ~ (^ 2AE ) 10^^'^ sec.~ 10 
years. More careful recent calculations that take into account the nuclear effects[13] are 
stated in terms of the following relation: 



'HVuc. — Rijnn) 



(5) 



where ~ 0.3 x 10^^ sec.~^, which yields tnuc. ~ 3 x 10^^ yrs. Present experimental limits 
on tnucI. from various nucleon decay searches e.g. at Kamiokande, Sudan, 1MB, Super- 
Kamiokande, SNO are given below in table I with the value of r„_^ noted there[16]. 



Expt 


source of neutrons 


TNuciiyrs) 


'^osc. (sec) 


Soudan 




0.72 X 10^2 


1.3 X 10^ 


Frejus 




0.65 X lO-'^^ 


1.2 X 10^ 


Kamiokande 


1%Q 


0.43 X 10^2 


1.2 X 10^ 


Super-K 


WQ 


1.77 X 10^2 


2.3 X 10^ 
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Cold Magnetic Vacuum 




L ~ 200 - 500 m 

Figure 1: Typical horizontal reactor set-up for n — n oscillation search 

Table Caption: Experimental lower limits on Tn-n and the nuclei used for deducing the 

limit. 

One can also anticipate a similar bound from the SNO experiment [20]. This implies 5m < 
10^^^ MeV. The nucleon decay searches for Tn-n become less efficient due to atmospheric 
neutrino induced background as we move to higher precision goals for Tn-n- It would therefore 
appear that in order to extend the Tn-n search much beyond the current limit, one needs to 
focus more on free neutron oscillation experiments rather than than nucleon decay search. 

We will see that present reactor neutron fluxes are precisely in the right range to probe 
these values of Tn-n that is compatibitible with matter stability as well as latest limits on 
proton decay. 



8 



2.2 Reactor Experiment to search for neutron oscillation 



The basic idea behind a reactor search for n — n oscillation is to have a cold neutron beam 
from a reactor (neutron speed f ~ 1 — 2 kilo meters/sec.) and have the beam pass through 
a high vacuum, demagnetized region for as long as feasible to a detector at the end of this 
"tunnel" . If one of the neutrons in the beam oscillates to an anti-neutron, this will annihilate 
in the detector and create a multi-pion signal with total energy of ~ 2 GeV. The need for 
demagnetization comes from the fact that the magnetic field of the Earth (~ 0.5 Gauss) 
will spht the neutron and anti- neutron energy levels by ~ 3 x 10~^^ GeV. If we take a 
typical flight time in the "tunnel" of about one sec, it will imply UnBtsimd x 10~^^ GeV 
sec which about a 1000 times h and hence does not satisfy the condition for free neutron 
oscillation. It is therefore necessary to reduce the magnetic field to the level of 10~^ — 10~^ 
Gauss. This is known technology that uses ^u-metal shielding. The typical horizontal set-up 
for such an experiment is given in Fig. 1 above. Under these conditions, one gets for the 
number of anti- neutrons for a given neutron beam the figure of merit of the experiment: 



where is the reactor flux; w„t= distance to detector; T running time. Maximum reactor 
fluxes for a 100 MW reactor is about ~ 10^^ — 10^^ neutrons per sec; for t — 0.1 sec. and 
T ~ 3 years can yield a hmit of 10^° sec. Using this technique, the first n — n search was 
carried out at the ILL laboratory in Grenoble, France and a lower limit of Tn-n > 8.6 x 10^ 
sec. was obtained[15]. New searches have recently been proposed which could improve this 
limit by two orders of magnitude [17]. Fig. 2 compares the free oscillation time to the nuclear 
instability life time (figure courtesy of Y. Kamyshkov Ref.[17]). It is clear that time of flight 
t anf the flux are two crucial factors that determine how precisely one can measure the 
n — n transition time. Clearly, the slower the neutrons the bigger is the time of flight. A 
possibility that is being pursued is to use ultracold neutrons with ~ 100 nano-eV kinetic 



^of events 




(6) 
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Figure 2: : 

Comparision of free neutron oscillation time vs nuclear instability life time 
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energies which move with speeds of order 4-5 m/sec and are trapped in a bottle[18] and may 
become competitive with cold reactor beam experiments in near future[19]. 



3 Operator analysis of AB = 2 processes and new physics 
scale probed 

Whenever a new process is used to probe physics beyond the standard model, a general 
question that is essential to know is the mass scale probed by it. A simple way to get 
a general idea about the scale probed by a physical process is to do an operator analysis 
i.e. assuming a particular symmetry and spectrum of the low energy theory, write effective 
higher dimensional operators for the process under consideration and see for what value of 
the mass scaling the operator, the process is observable. This argument has been a useful 
tool to probe physics at short distance scales probed by processes such as proton decay, 
neutrino mass as well as corrections to standard model observables. 

The method however has its limitations. For instance, often in these discussions, one 
uses the SM fermion spectrum but if there is a SM non-singlet new particle with a 100 GeV 
mass not discovered yet, the naive scale arguments can be misleading since new operators 
may appear. Similarly, the presence of unknown higher symmetries could also invalidate 
these arguments. We will apply this discussion to the AB — 2 operators below to see the 
scale probed by neutron-anti-neutron oscillation. 

3.1 Standard model operator 

In the standard model, the effective operator responsible for AB — 2 and AL — processes 
is given by: 

Oab=2 = ^[gQgg(d'^d')* + iiWiiW] (?) 
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The strength of this operator Gab=2 — jg?- The n — n mixing mass can be deduced from 
this operator by simple dimensional analysis to be: 

6mn-n ^ cGaB=2A-qcD = ^^jf (S) 

where we chooze c to be of order one. Attempts have been made in the past to estimate c 
using bag as well as other phenomenological models for hadrons[21]. Taking the best current 
lower bound on Tn^n from ILL reactor experiment [15] which is 10^ sec. and comparable 
bounds from nucleon decay search experiments [16], one can then obtain a lower limit on M 
to be around 10-300 TeV depending on other couplings in a theory. If proposals to improve 
the precision of this search by at least two orders of magnitude [17] are carried out, , then 
in the context of the standard model particles, the mass scale probed will be around a 1000 
TeV. Thus we see that neutron-anti-neutron oscillation will probe physics at scales far below 
the GUT scales that proton decay will probe. An important question to ask is whether the 
scale probed can be higher in the presence of new physics around the TeV scale. We discuss 
this below. 

3.2 Supersymmetry and enhancement of AB = 2 operator 

In the presence of supersymmetry at the TeV scale, there are new particle such as squarks 
and sleptons with TeV scale mass. They can then enter the effective AB = 2 operator. This 
kind of effect is familiar from discussions of supersymmetric grand unified theories where the 
presence of super-partners with TeV scale masses can drastically alter the operator analysis 
for nucleon decay. For instance the dimension six operator responsible for proton decay in 
non-susy models goes like QQQL/M'^ whereas the presence of squarks with TeV or sub- TeV 
masses, the dominant operator becomes QLQQ/M, which has reduced power dependence 
on mass and is well known to put severe constraints on grand unified theories. 
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A typical leading operator for n — n oscillation for supersymmetric case is: 

Oab=2 = ■^u'^d'^d-u-d'^d'^ (9) 

Note that the power dependence on the seesaw scale (or B-L breaking scale) has now consid- 
erably softened. The conversion of susy particles to SM particles brings some suppression; 
but the the overall impact is that simple power counting arguments change. They imply 
that if this is the leading operator one can probe the seesaw scale upto 10^ GeV with n — n 
oscillation times of 10^° — 10^^ sec. The power dependence in fact softens even further if 
there are new color sextet particles at or below the TeV scale. For instance, if there is a 
scalar diquark coupled sextet field of A^c^c type, the leading operator becomes: 

Oab=2 = ^«A„c„cJ-J'^ (10) 

and the scale reach oin — n goes up to 10^^ GeV. If on the other hand there is a field A^cdc, 
the leading order operator becomes: 

Oas=2 = -^d^d^Auc^cAucdc (11) 

increasing the scale reach to the GUT scale. The question then arises whether there are 
plausible models where this can happen. Below I give an example of models which have TeV 
scale Auc„c fields naturally so that an 1/M^ dependence on the seesaw scale is quite natural. 
Prom this, we see that scale of observable n — n physics can be anywhere from 300 TeV to 
10^3 TeV 

4 Neutrino mass physics probed hy n — n oscillation 

To see what new physics is revealed by the search for n — n oscillation, we give several 
examples where low scale for seesaw manifests in the non-leptonic sector of the theory as 
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n — n oscillation. This happens due to the connection of seesaw to local or global B-L 
symmetry. There could also be other kind of beyond the standard model scenarioss where 
extra space dimensions curled up to TeV~^ length scales and where additional dynamical 
inputs can lead to observable n — n oscillation. In the former case, it will throw important 
light on the physics associated with neutrino mass and we discuss this below. 

There are various reasons to think that B-L symmetry is an intimate feature of physics 
beyond the standard model. Some of them are as follows: 

• (i) The seesaw mechanism for understanding small neutrino masses requires the intro- 
duction of right handed neutrinos [5]. In the presence of three right handed neutrinos it 
is more natural for B-L symmetry which was a global symmetry of the SM Lagrangian 
to become a local symmetry; 

• (ii) An inherent aspect of seesaw is the Majorana mass of the right handed neutrino 
that breaks the B-L symmetry providing one way to understand why the seesaw scale 
is so much less than the Planck scale; 

• (iii) A third reason appears once one admits the presence of supersymmetry at the 
TeV scale, as is widely believed and requires the hghtest SUSY particle (LSP) to be 
naturally stable in order to be the dark matter candidate of the universe. The simplest 
way to have SUSY LSP naturally stable is to have B-L as a symmetry of physics beyond 
MSSM[7]. 

If indeed B-L symmetry is present in nature, there are several questions that immediately 
come to mind: 

• (a) Is it a global or local symmetry ? 

• (b) is it a broken or unbroken symmetry ? 
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• (c) if it is broken, what is the breaking scale and what is the associated physics ? 



I will now argue that search for the process of ncutron-anti-neutron oscillation [8, 9, 10] can 
provide a partial answers to some of these important questions: in particular the question of 
scale of B-L breaking and associated expanded gauge symmetry and hence physics associated 
with the origin of neutrino mass. The point is that since a Ai? = 2 operator breaks B-L 
by two units, it is natural to associate the scale associated with the Ai? — 2 operator also 
with the scale of B-L breaking since neither B nor L are separate symmetries of the standard 
model but rather the symmetry is B—L . Since seesaw mechanism also breaks B-L symmetry, 
M could also be the seesaw scale. In fact, in the context of models such as those based on 
SU{2)l X SU{2)r X 5'C/(4)c[14], the same mass scale M is responsible for both processes as 
was first shown in Ref.[10]. Another purely group theoretical way to see this is to note that 
in left-right symmetric models, there is a relation between the electric charge and the B-L 
generator as follows: 

Q ^ hL + ISR + (12) 

For distance scales shorter than the electroweak scale, we have AI^l — and electric charge 
conservation then implies that AI^r = — implying that parity violation (or AI^r = 1) 
implies not only that AL — 2 i.e. seesaw mechanism for neutrino mass but in the non- 
leptonic sector it implies AB — 2 which is of course n — n oscillation. 

Thus the search for n — n oscillation could not only illuminate the nature of the important 
symmetry, B-L but could also be one way to unravel the mystery of the seesaw mechanism 
that is expected to be major player in the physics of neutrino mass as well as the origin of 
parity violation in Nature. 
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5 SU{2)l X SU{2)r X SU{A)c model with light diquarks 



In this section we demonstrate by an explicit example how neutrino mass and n — n are 
intimately connected in this partial unification theory. We first recapitulate some elementary 
facts about the model. 

The quarks and leptons in this model are unified and transform as ip : (2, 1, 4) © ■^'^ : 
(1, 2, 4) representations of SU (2)^ x SU (2)^^ x SU (4)c. For the Higgs sector, we choose, 0i : 
(2, 2, 1) and ^15 : (2, 2, 15) to give mass to the fermions. The A'^ : (1, 3, 10) ©A^ : (1, 3, 10) 
to break the B — L symmetry. There are color sextet scalar diquark fields contained in the 
A"^ : (1, 3, 10) multiplet, which will play a crucial role in generating n — n oscillation. We 
first discuss the non-Osupersymmetric version of the model in next subsection and follow it 
up by the supersymmetric generalization where the scale reach of n — n goes up all the way 
to 10^° GeV or so. 

5.1 Non-supersymmetric version 

The Higgs potential for this model can be written as : 

V{^oAl^ Ar,l) = K) + Kn (13) 

where 

Vm = A^e'''^'^^e^'''^''^''''AR,^^,A,,,A,,,A,,, + {R^L)+ h.c. (14) 

This term is crucial for generating baryon number non-conservation in this theory and with- 
out it, the model cannot explain the origin of matter in the universe (as we see below). We 
will see that the same term is responsible for the AB — 2 n — n oscillation. 

In order to discuss neutrino masses and n — n oscillation, let us write down the Yukawa 
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couplings in the model: 



Ly = + /ll5*L^15*ii + f^\tC-^^R^R +(R^L)+ H.C. (15) 

Using standard spontaneous breaking via the vev of (1, 3, 10) to the standard model gauge 
group and $1,15 to give mass to both to charged fermions and Dirac mass to neutrinos, 
we implement the seesaw mechanism for neutrinos. It was shown in ref.[10] that the same 
(1, 3, 10) via the potential term leads to n — n oscillation via the diagram in Fig. 3. The 
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Figure 3: A Feynman diagram contributing to n — n oscillation as discussed in Ref. 8. 
Another class of graphs involves two A^c^^c and one A^c^^c scalar exchanges. 



strength of n — n oscillation is given by: 

Gab=2 — 



Ml 



(16) 



We see the fifth power dependence on the seesaw scale since the diquarks are expected to be 
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at seesaw scale by the usual naturalness considerations. Also note that even after symmetry 
breaking this model has a Z2 symmetry given by (—1)^. 

It has recently been shown[22] that this model also explains the baryon asymmetry of 
the universe via the mechanism of post-sphaleron baryogenesis to be discussed below[23]. 

5.2 Supersymmetric version 

We now discuss the supersymmetric version of this model where the triplet fields doubled to 
cancel the gauge anomalies. We also add a, B — L neutral triplet Q : (1, 3, 1) which helps to 
reduce the accidental global symmetry of the model and hence the number of light diquark 
states. The superpotential of this model is given by: 

W ^ Wy + Whi+ Wh2+ Wh3 (17) 

where 

Whi = XiS{A'A'-M^) +XcA'A'fl + //^Tr (0,0^) (18) 



(A-A-)^ , , (A-A-)(A-A-) 



. (.0, 

(21) 



Wy = h2P<P,r + hi5^<Pi5r + fr^'r- (22) 

Note that since we do not have parity symmetry in the model, the Yukawa couplings hi and 
hi5 are not symmetric matrices. When Ab = 0, this superpotential has an accidental global 
symmetry much larger than the gauge group[25]; as a result, vacuum breaking of the B — L 
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symmetry leads to the existence of light diquark states that mediate N N oscillation and 
enhance the amplitude. In fact it was shown that for (A'^) ~ (A'^) 7^ and (Q) 7^ and all 
VEVs in the range of 10^^ — 10^^ GeV, the light states are those with quantum numbers: 
Ayc^c. The symmetry argument behind is that [25] for As = 0, the above superpotential 
is invariant under f/(10,c) x SU{2,c) symmetry which breaks down to U{9,c) x U{1) when 
(A^cj^c) = vbl 7^ 0. This results in 21 complex massless states; on the other hand these 
vevs also breaks the gauge symmetry down from 5'C/(2)r x 5'C/(4)c to SU{S)c x U{1)y- This 
allows nine of the above states to pick up masses of order qvbl leaving 12 massless complex 
states which are the six A^c^c plus six A^c„c states. Once 7^ and is of order 10~^ — 10~^, 
they pick up mass (call Muc„c) of order of the elctroweak scale. 

5.3 A new diagram for neutron- anti- neutron oscillation 

To discuss n ^ n oscillation, we introduce a new term in the superpotential of the form[10]: 

Wab=2 = ^e'^''^'^''^'e^^^'^A^^,A^,,A^,,A^,,, (23) 

where the /i, z/ etc stand for SU{4:)c indices and we have suppressed the SU{2)r indices. 
Apriori could be of order Mpf, however the terms in Eq.(2) are different from those in 
Eq. (4); so they could arise from different a high scale theory. The mass is therefore a 
free parameter that we choose to be much less than the Mpe. This term docs not affect the 
masses of the Higgs fields. When A^c,^c acquires a VEV, AB — 2 interaction are induced 
from this superpotential, and n n oscillation are generated by two diagrams given in 
Fig. 3 and 4. The first diagram (Fig. 3) in which only diquark Higgs fields are involved 
was already discussed in [10] and goes like Gn^n ~ m^^^^m'-^^m ' Taking M^c^c ~ 350 GeV, 
^d'^d'^ ~ A'l'sL and Ma ~ vbl as in the argument [25], we see that this diagram scales 
like Vb'Iv~1. In ref.[24] a new diagram (Fig. 4) was pointed out which owes its origin to 
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supersymmetry. We get for its contribution to Gab=2- 

fflVBL 



Gr. 



9l 



(24) 



167r2 M2e„cM|e,cMsusYM, " 

Using the same arguments as above, we find that this diagram scales hke v^1v~l which is 
therefore a significant enhancement over diagram in Fig.3 when vbl ^ Vwk- In order to 
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Figure 4: The new Feynman diagram for n — n oscillation. 



estimate the rate for n ^ n oscillation, we need not only the different mass values for which 
we now have an order of magnitude, we also need the Yukawa coupling /n. Now /n is 
a small number since its value is associated with the lightest right-handed neutrino mass. 
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However, in the calculation we need its value in the basis where quark masses are diagonal. 
We note that the n — n diagrams involve only the right-handed quarks, the rotation matrix 
need not be the CKM matrix. The right-handed rotations need to be large e.g. in order to 
involve /33 (which is 0(1)), we need {V^'%i to be large, where vi''''^^W^^dV^''^'^ = l^^^". 
The left-handed rotation matrices V^"'''-' contribute to the CKM matrix, but right-handed 
rotation matrices V^''^^ are unphysical in the standard model. In this model, however, we 
get to see its contribution since we have a left-right gauge symmetry. Let us now estimate 
the time of oscillation. When we start with a /-diagonal basis (call the diagonal matrix 
/), the Majorana coupling /n in the diagonal basis of up- and down-type quark matrices 
can be written as (VJ/Vr)!! ~ {V^iThs- Now fss is 0(1) and can be ~ 0.6, so fu is 
about 0.4 in the diagonal basis of the quark matrices. We use Msusy, MuCu<= ~ 350 GeV and 
vbl ~ 10^^ GeV. The mass of A^cdc i.e. M^c^^c is 10^ GeV which is obtained from the VEV 
of Q : (1, 3, 1). We choose M* ~ 10^^ GeV. Putting all the above the numbers together, we 
get 

Gn^n ^ 1 • 10-'° GeV-^ (25) 

Taking into account the hadronic matrix element effect, the n — n oscillation time is found to 
be about 2.5 x 10^° sec which is within the reach of possible next generation measurements. 
If we chose, ~ Mp^, we will get for Tn-n ~ 10^^ sec. unless we choose the M^d^c to 
be lower (say 10^ GeV). This is a considerable enhancement over the nonsupcr symmetric 
model of [10] with seesaw scale of 10^^ GeV. We also note that as noted in [10] the model is 
invariant under the hidden discrete symmetry under which a field X — > e^'^^^X, where Bx 
is the baryon number of the field X. As a result, proton is absolutely stable in the model. 
Furthermore, since R-parity is an automatic symmetry of MSSM, this model has a naturally 
stable dark matter. 
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6 Baryogenesis and n — n oscillation 



In the early 1980's when the idea of neutron-anti-neutron oscillation was first proposed in 
the context of unified gauge theories, it was thought that the high dimensionality of the 
AB 7^ operator would pose a major difficulty in understanding the origin of matter in 
the Universe. The main reason for this assessment is that the higher dimensional operators 
remain in thermal equilibrium until late in the evolution of the universe without contradicting 
any low energy observations. This is because the thermal decoupling temperature T* for 
such interactions goes roughly like vbl ( if^ ) which can be in the range of temperatures 
where B-l-L violating sphaleron transitions are in full thermal equilibrium, say for example 
Vbl ^ 10 — 100 TeV, as required for the case of observable n — n oscillation. They will 
therefore erase any baryon asymmetry generated in the very early moments of the universe 
(say close to the GUT time of 10~^^ sec. or so) in then prevalent baryogenesis models. Even 
though GUT baryogenesis models are no more popular, the same argument will apply to 
any other high temperature baryogenesis mechanism. 

In models with observable n — n oscillation therefore, one has to search for new mech- 
anisms for generating baryons below the weak scale. In this section, we discuss such a 
possibility [2 3] which was discussed in a recent paper. As we see below, it is ideally suited for 
embedding into the G224 model discussed in Ref.[10] and leads to an interesting side "bonus" 
that it puts an upper limit on the neutron-anti-neutron oscillation time Tn-n. 

As an illustration of the way the new mechanism operates, let us assume that there is a 
complex scalar field that couples to the AB — 2 operator as follows i.e. Lj — Su'^dPd'^u'^d'^d'^ / , 
where the scalar boson has mass of order of the weak scale and B = 2. When < S >^ 0, this 
interaction leads to baryon number violation by two units and observable n — n transition 
if M is in the few hundred to 1000 GeV range. Writing S = -j^{vbl -\- Sr -\- Sj), we see 
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that the direct decay of Sr involves both AB — ±2 final states and thereby satisfies the first 
requirement for baryogenesis. 

The first point to note is that the high dimension of Lj allows the scalar AB ^ decay 
to go out of equilibrium at weak scale temperatures. This satisfies the out of equilibrium 
condition given by Sakharov for generating matter-anti-matter asymmetry. To see this, 
let us give some examples: if we chose the proton decay operators such as QQQL, the 
decoupling temperature consistent with present experimental bounds on proton life time 
would be around 10^^ GeV or so. So to apply our mechanism, we need to consider higher 
dimensional operators. A typical example of such an operator that we will focus on in this 
paper is of type: u'^(FdFu'^dF(F/M^. As note earher, this operator leads to n — n oscillation 
and for this to be in the observable range, the associated mass scale has to be in the few 
TeV range. The decupling temperature can then be easily estimated from the formula 
Trf ~ M {jjj^^ ^ ~ 10~^'^M. For appropriate values of M in the range of interest can be 
below 100 GeV so that the sphalerons have gone out of equilibrium and baryogenesis follows. 
To make these ideas more concrete, below we give an explicit example[23]. 

We consider an effective sub- TeV scale model that gives rise to the above higher dimen- 
sional operator for n <-> ri oscillation. It consists of the following color sextets, SU (2) l sin- 
glet scalar bosons {X, Y, Z) with hypercharge — |, +|, +| respectively that couple to quarks. 
These states emerge naturally if this low scale theory is embedded into the SU (2) l x SU (2)^ x 
SU{A)c (denoted by G224) model described in the previous section. The scalar field S can be 
identified with the complex field A^c^c whose vev breaks B-L symmetry by two units as is the 
case for us and Sr = i?e(A°Ci^c[22]. The X, F, Z fields can be identified with the color sextet 
partners of Ai,c,^c in the (1,3,10) multiplet i.e. the A^cj^c, Auc^c and Auc^c respectively. In 
the X, Y, Z notation, one can write down the following standard model invariant interaction 
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Lagrangian: 

Lj = hi^Xd1d] + fi^Yu1u] + (26) 
gijZ{u1d] + u''^df)+ XiSX^Y + MSXZ^ 

The complex scalar field 5" clearly has B — L = 2. To see the constraints on the parameters 
of the theory, we note that the present limits on Tn-n > 10*^ sec. implies that the strength 
Gn-n of the the AB — 2 transition is < 10~^^ GeV~^. Prom the above figure, we conclude 
that 

For Ai_2 ~ /ill ~ /ii ~ S'li ~ 10~^, we require Mi ~ Mx,y,z — 1 TeV to satisfy this 
experimental bound. In our discussion of generic models of this type, we will stay close 
to this range of parameters and see how one can understand the baryon asymmetry of the 
universe. The singlet field will play a key role in the generation of baryon asymmetry. We 
assume that < 5" Mx and Ms^ ~ 100 — 500 GeV. It can then decay into final states 
with B = ±2 i.e. six quarks and six anti-quarks. 

Note that the Lagrangian of Eq. (26), leads to tree level contribution to flavor changing 
neutral current processes such as K — K, B — B mixings etc and that will restrict the form of 
the flavor structure of the coupling matrices suitably. These constraints can also be satisfied 
in a realistic 6*224 embedded model which also explain the neutrino masses and mixing[22]. 
This is an important consideration since the magnitude of the baryon asymmetry depends 
on this detailed flavor structure. 

On the way to calculating the baryon asymmetry, let us first discuss the out of equilibrium 
condition. As the temperature of the universe falls below the masses of the X, Y, Z particles, 
the annihilation processes XX d'^d'^ (and analogous processes for Y and Z) remain in 
equilibrium. As a result, the number density of X, Y, Z particles gets depleted as e~^^^'^ 
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and only the Sr particle survives along with the usual standard model particles. One of the 

primary generic decay modes of Sr is Sr u'^d^d'^u^d'^d^. There could be other decay modes 
which depend on the details of the model and one has to ensure that the rates to other 
modes are smaller than the six quark mode. 

The generic chain of events leading to baryogenesis in this model is the following. At 
T ~ Ms^, the decay rate is smaller than the Hubble expansion rate. As the Universe cools 
below this temperature, the decay rate remains constant whereas the expansion rate of the 
universe is slowing down. So at a temperature far below Ms, S will start to decay when 
the decay rate F ~ with is given by: 

Since the corresponding epoch must be above that the QCD phase transition temperature, 
this puts a constraint on the parameters of the model. Typically, we need to have Ms / Mx ~ 
0.5 or so due to high power dependence of on this ratio. This implies that the X, Y, Z 
masses cannot be arbitrarily high, since the heavier these particles are, the lower will be. 
We expect this upper hmit to be in the TeV range at most as we show below . 

It is well known that baryon asymmetry can arise only via the interference of a tree 
diagram with a one loop diagram. The tree diagram is clearly the one where 5" — > 6g. There 
are however two classes of loop diagrams that can contribute: one where the loop involves the 
same fields X, Y and Z. A second one involves W-exchange, which involves only standard 
model physics at this scale (Fig. 6). We find that the second contribution can actually 
dominate. In fact, in the G224 embedding of the model, the first diagram vanishes. The 
second diagram also has the advantage that it involves less number of arbitrary parameters 
and the source of CP violation in this case is the same as the CKM CP violation present in 
the standard model. The observed baryon asymmetry is related to the primordial baryon 
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asymmetry cb as r)B — (^pI^^^cb) where 



We find that 



en ^ 



nsr{S^6q)^T{S^6q) 



r{S) 



(29) 





Figure 5: Tree level Feynman diagram for S decays. 
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Figure 6: One loop diagram for S decays. 



, vertex 



a2 lmTr[gg^M^VMag^M^V*Md] 



(30) 



Br 4 Tr{g^g)M^Ml 

The magnitude of the asymmetry depends on Ms as well as the detailed profile of the various 
coupling matrices h,g, f and we can easily get the desired value of the baryon asymmetry by 
appropriately choosing them. An important consequence of the above equation regardless 
of the details is that if Ms is much bigger than about 500 GeV, baryon asymmetry becomes 
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very small. This implies that Mx,y,z in turn should not be much larger than a TeV, implying 
that they can be accessible at the LHC This also implies that the strength oln — n transition 
has a lower limit or an upper limit on Tn^n. 

There is a dilution of the baryon asymmetry arising from the fact that -C Ms since 
the decay of S also releases entropy into the universe. Thius dilution factor is 

^ (0.32g4Ta)Taf/' 

" (O.I2M5 + 0.32aTi3)3/4 ^""-"^ 
Since the decay rate of the S boson depends inversely as a high power of Mxy, higher X, Y 
bosons would imply that ~ i7 is satisfied at a lower temperature and hence give a lower 
ds- For our preferred choice of parameters i.e. Mx ~ 2M5 ~ TeV, we find ds ~ 20%. Also 
for choice of the coupling parameters Xr^fr^h^g^^ 10~^, and Ms — 200 GeV we find 
Tn-n < 10^° sec. 

We also note that if the model is embedded into a G224 group, all three coupling matrices 
h, g, f become equal; there are also some constraints on the gauge boson spectrum in this 
model[22]. As noted above, an interesting consequence of adequate baryogenesis is that there 
are new sub-TeV- TeV scale color sextet scalar bosons which can be observed at LHC [26]. 

7 Extra dimensional models with observable n — n os- 
cillation 

The possibility that there may be extra compact space dimensions has been under extensive 
theoretical as well as experimental investigation inspired by the fact that superstring the- 
ories do predict the existence of such models. In a large class of interesting models, these 
dimensions are believed to be of order TeV~^ size[27, 28] making their manifestations in 
many physical situations of experimental interest not only for colliders but also for low en- 
ergy experiments. One class of such possibilities includes the breakdown of baryon number 
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by two units. We note two such discussions: 

(i) In ref.[29], it was noted that if there are miUimeter size extra dimensions, one way to 
understand the observed fermion mass hierarchies would be to distribute the standard model 
fermions in the bulk. In such a scenario, their Yukawa couplings would be suppressed if they 
are further apart [30]. Thus in the context of a 5- or 6 dimensional scenario, observed mass 
hierarchies pretty much fix the locations of the various standard model fermions in the bulk. 
One can then ask the question as to what happens to higher dimensional operators such as 
those giving proton decay and n — n oscillation. Since these are low scale models, one has to 
separate the quarks and leptons sufficiently to suppress these operators since their strength is 
given by the overlap integral of the various wave functions in the extra dimension. The n — n 
operator however involves only the right handed quarks and their locations is generally fixed 
by the mass hierarchy considerations. So there is no freedom to suppress their strengths. 
The strength of such operators was calculated in the context of six dimensional models in 
Ref.[29] and found to be in the observable range for extra dimension sizes of order 45 — 100 
TeV, which are of the same order as being contemplated for phenomenological viable models 
of this kind. These considerations could also be applied to the case Randall-Sundrum models 
where the SM fermions are also generally expected to be distributed in the bulk. 

(ii) A second class of models were discussed in ref.[31], where general considerations of 
global charge separation from the branes is used to argue that one cannot have proton decay 
but observable n — n oscillations due to the fact that neutrons are both color and electric 
charge neutral. Such baryon number nonconservation occurs due to brane fluctuations which 
can create baby branes that carry baryon number if they are color and charge neutral. 
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Figure 7: neutron oscillation induced in R-parity breaking MSSM 

8 R-parity breaking models 

Finally, there have been explorations of n — n oscilation to supersymmetry. As is well 
known, in minimal supersymmetric extensions of the standard model (MSSM), a baryon 
number violating operator of the form u'^d'^d'^ is allowed in the supcrpotential. By itself it 
can lead to Ai? = 2 operator of the form {u'^d^d'^y at tree level via gluino exchange (Figure 
7). Typical strength of such operators arising from this new interaction is given[32] by: 

^^^=^ - ~MIM~ ^^^^ 
Note that due to color anti-symmetry, at the tree level, this interaction leads A — A oscillation 
rather than n — n oscillation. However once radiative corrections are included, this operator 
via flavor changing effects gives rise to neutron oscillation [33]. This puts constraints on 
the R-parity violating couplings of the A" type. It must be noted that if these couplings 
are present in combination with R-P violating operators of type QLd'^, then it leads to 
catastrophic proton decay rates[34]; so in discussing AB = 2 transitions, we are assuming 
that QLd'^ operators are absent. 



9 Conclusion 

In summary, there is now a widely held belief that neutrino mass provides a clue to the 
existence of a new local B-L symmetry of nature. A key question in particle theory beyond 
the standard model is then to find out the scale and detailed nature of physics associated 
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with this new symmetry. The situation is similar to the late 60's when the hints in favor 
of the existence of a fundamental local SU (2) symmetry of weak interaction were becoming 
compelling that eventually led to the triumphant gauge revolution in the 1970's. They were 
subsequently confirmed by the discovery of neutral currents, followed by the discovery of the 
W and Z boson. The two paths before us right now arc (i) whether the B-L forces responsible 
for neutrino mass represent physics close to the gravity scale of 10^® GeV as exemplified, say 
by the grand unified theories or (ii) scales multi-TeV range far below the GUT scale with 
signals visible to the Large Hadron Colliders. One direct way to explore this new physics is 
to to carry out an experimental search for n — ri oscillation. As noted the discovery oin — n 
oscillation will have far reaching implications not only for the nature of forces and matter 
but will completely alter our thinking about such cosmological issues as the origin of matter, 
nature of dark matter etc. 

The main points stressed in this review are that: (a)there are a wide class of models 
e.g. SU{2)l X SU{2)ji x SU{4:)c with type II seesaw, models based on TeV scale extra 
dimension models where n — n oscillation is in a range accessible to currently planned reactor 
experiments and they will extend the existing limit on the strength of this process to an 
exciting range; (b) Models that predict observable n — n oscillation provide a new way to 
understand the origin of matter where matter creation moment is "much" later than currently 
contemplated models such as leptogenesis or GUT scale baryogenesis. More importantly, 
unlike the other mechanisms for baryogenesis, this mechanism can be tested in current 
collider experiments as well as neutrino experiments in some specific realizations. 

On the experimental front, searches for n — n transition using the proton decay type 
set-ups are back ground limited [35] and cannot be used to push the limits for Tn-n too much 
beyond the present limit. Also an interesting class of theories that gives rise to n — n induced 
baryogenesis lead to new coloe sextet particles observable at LHC[26]. 
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Finally, it is also worth noting that there is a related phenomenon involving neutrons 
i.e. neutron-mirror neutron oscillation that can also be searched for in experimental set- 
ups similar to that searching for n — n oscillation [36]. This new class of oscillations have a 
much weaker experimental limit on its strength and could be improved by longer "baseline" 
searches being contemplated for n — n oscillation. 

I would like to thank K. S. Babu for many discussions on theoretical aspects of the topic 
and Y. Kamyshkov, G. Greene, Mike Snow and Albert Young for many discussions on the 
experimental prospects forn — n oscillation. I am grateful to K. S. Babu and Y. Kamyshkov 
for carefully reading the manuscript and suggesting improvements. This work is supported 
by the National Science Foundation grant no. Phy-0652363. 
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